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Abstract — We study interference channels (IFC) where inter- 
action among sources and destinations is enabled, e.g., both 
sources and destinations can talk to each other. The interaction 
can come in two ways: 1) for half-duplex radios, destinations 
can talk back to sources using either simultaneous out-of-band 
(white spaces) transmission or in-band half -duplex transmission; 
2) for full-duplex radios, both sources and destinations can 
transmit and listen in the same channel simultaneously. The 
flexibility afforded by interaction among sources and destinations 
allows for the derivation of interference alignment (IA) strategies 
that have desirable "engineering properties": insensitivity to the 
rationality or irrationality of channel parameters, small block 
lengths and finite SNR operations. We show that for several 
classes of interference channels the interactive interference align- 
ment scheme can achieve the optimal degrees of freedom. We 
also give simulation results on the finite SNR performance of 
interactive interactive alignment for 3 -user and 4-user half- 
duplex interference channels. 

I. Introduction 

Due to the broadcast nature of wireless communication, 
an active transmitter will create interference to all other 
receivers using the same channel (in time and frequency). 
This is a basic challenge in supporting concurrent reliable 
wireless communications. Traditional methods have involved 
power control, scheduling and antenna beamforming/nulling 
to manage interference. Interference alignment is a novel 
interference management technique which makes concurrent 
wireless communications feasible 0], 0. 

Since its introduction, substantial research has been con- 
ducted in understanding the gains of interference alignment. 
From a theoretical perspective, the focus has been mainly 
in understanding the degrees of freedom (first term in the 
high SNR approximation of the information theoretic capacity 
region). A seminal result, shown in (T), is that we can 
communicate at a total ~ degrees of freedom for nearly all 
time-varying/frequency selective i^-user interference channel 
(this is also an upper bound). To get to § degrees of freedom, 
the scheme employed (vector space interference alignment) 
requires the channel diversity to be unbounded; in fact, one 
needs the channel diversity to grow like K 2K 0. The 
diversity order required is huge even for modest number of 
users K (such as 3 and 4). This is definitely so in a practical 
communication system, where there is hardly enough channel 
diversity to make such a vector space interference alignment 
scheme feasible. 

A practically relevant problem is to understand the funda- 
mental degrees of freedom for a fixed deterministic channel. 
For fully connected channel matrix H, the total degrees of 



freedom are upper bounded by ~ and J4) shows that this 
upper bound is achievable for almost all channel matrices H 
using a coding scheme based on Diophantine approximation. 
However, this result is limited in two ways. First, the coding 
scheme is very sensitive to whether entries of H are rational 
or irrational. Second, although it is provable y degrees of 
freedom are achievable for almost all H, for a given H, in 
general it is not known what are the optimal achievable degrees 
of freedom. shows that the y result for almost all H 
can be derived using Renyi information dimension. Again, 
the result is sensitive to whether entries of H are rational 
or irrational, and for fixed channel matrix H, in general 
the optimal achievable degrees of freedom is unknown. The 
recent works of @, Q address this issue to a good extent 
for the case of the two-user X channel and the symmetric 
K-nsei interference channel respectively, but the engineering 
implication of the proposed coding schemes remains unclear. 

Therefore, despite significant theoretical progress on the 
if -user interference channel problem, it is still unclear how 
to make interference alignment practical. The drawbacks of 
existing schemes may be inherent to the channel model which 
assumes sources can only transmit and destinations can only 
listen, while in practice radios can both transmit and re- 
ceive. We study new channel models where interaction among 
sources and destinations is enabled, e.g., both source and 
destination can talk to each other. The interaction can come 
in two ways: 1) for half-duplex radios, destinations can talk 
back to sources using simultaneous out-of-band transmission 
or in-band half-duplex transmission; 2) for full-duplex radios, 
both sources and destinations can transmit and listen in the 
same channel simultaneously (8). Although shows that 
for interference channel, relays, feedback, and full-duplex 
operation cannot improve the degrees of freedom beyond 
~, we demonstrate that the interaction among sources and 
destinations enables flexibility in designing simple interference 
alignment scheme and in several cases achieves the optimal 
degrees of freedom. 

Our main contribution is to propose a simple interfer- 
ence alignment scheme by exploiting the interactions among 
sources and destinations, and prove that the scheme can 
achieve the optimal degrees of freedom for several classes of 
interference channels, including half-duplex 3-user IFC, full- 
duplex 3-user and 4-user IFC, and full-duplex 4-user MIMO 
IFC. One specific aspect of our model, namely, feedback using 
the reciprocal interference channel, has been considered in 
prior work iflOll . In this work, we improve on this state-of-the- 



art in two ways: we prove new results for this specific model 
and also generalize this model to exploit more general modes 
of interaction in order to achieve interference alignment. In 
addition, we do extensive numeric simulations and show the 
proposed interactive interactive alignment scheme also works 
for some other classes of IFC empirically. We use tools 
from algebraic geometry to show why success of numeric 
simulations can suggest that the scheme should work well for 
almost all channel parameters in a rigorous way. 

The paper is organized as follows. We describe the channel 
model in Section and give some motivating examples in 
Section III We present the interactive communication scheme 
for half-duplex radios and interference alignment conditions 



Sources 



in Section IV and our main technical results from algebraic 
geometry on the interference alignment feasibility in Section 
|V| Section VI and Section VII study interactive interference 
alignment for half-duplex and full-duplex K-usei interference 
channels with small K, respectively. Simulation results on the 
finite SNR performance evaluation are presented in Section 



VIII Section IX discusses a general multi-phase interactive 



communication scheme for i^-user interference channel with 
large K, and lists several open problems. Section [X] concludes 
this paper. 

II. System Model 

Consider a K-user interference channel with 2K ra- 
dios, where radios Sx, S2, . . . , Sjc are sources and radios 
t\,t2, ■ ■ ■ ,tx are destinations. For each 1 < i < K, source Si 
wants to send an independent message Wi to destination ti. 

Let H £ C denote the forward channel matrix from 
sources to destinations, and the input and output signals of the 
forward channel are related as 



K 



] = y^-HjjXj[n] + Zi[n], 



(1) 



where n is the time index, yi[n] is the signal received by 
destination ti at time n, Xj [n] is the signal sent out by source 
Sj at time n, Zi [n] is the channel noise, and ify is the channel 
coefficient from source sj to destination tj. 

The above is a canonical channel model for K-nsei inter- 
ference channel. 

In this work, we consider two channel models where inter- 
action among sources and destinations can be enabled. 

A. Interaction from destinations to sources via reverse trans- 
mission for half-duplex radios 

Since in practice radios can both actively talk (send signals) 
or passively listen (receive signals), we assume destinations 
can also actively talk back to sources either via simultaneous 
out-of-band (white spaces) transmission or via in-band half- 
duplex transmission. Let G £ C KxK denote the feedback 
channel (or reverse channel) matrix from destinations to 
sources. Similarly, the input-output relation of the feedback 




Fig. 1. System Model For Half-duplex Radios: s\, S2, ■ ■ ■ , sjj- are sources 
and ti, t2, ■ ■ ■ , tn are destinations. H is the forward channel matrix from 
sources to destinations, and G is the feedback channel matrix from destina- 
tions to sources. 



channel is 



fi[n] 



K 

E 



GijVjln] + Zi[n], 



(2) 



where n is the time index, Vj[n] is the signal sent out by 
destination tj at time n, fi [n] is the signal received by source 
Si at time n, Zi[n] is the channel noise, and Gij is channel 
coefficient from destination tj to source Sj. 

If the destinations use in-band half-duplex transmissions, 
due to the reciprocity of wireless channels, we have 



G = H T , 

where H T denotes the transpose of H. 



(3) 



B. Interaction via transmitting and receiving simultaneously 
for full-duplex radios 

Full-duplex radios can both send and receive signals using 
the same channel simultaneously, and this capability naturally 
enables the interaction among all radios in the network. 

We assume that all sources and destinations have full-duplex 
antennas, and all nodes can transmit and receive signals using 
the same channel simultaneously. Let H £ C KxK be the 
channel matrix from sources to destinations, U £ <C KxK be 
the channel matrix among sources, and W £ C KxK be the 
channel matrix among destinations. The input-output relation 
of this full-duplex communication channel is 



K 



K 



K K 

M n ] = E h m v A ti \ + E u a x A n \ + ^Wi ( 5 ) 

j'=i i=i 

where n is the time index, xj [n] and vj [n] are the signals sent 
out by source Sj and destination tj at time n, respectively, 
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• Phase 1: All sources send their independent symbols 
Xi,X2,X3 simultaneously. Then destinations receive 



Fig. 2. System Model For Full-duplex Radios: si, S2, ■ ■ ■ , are sources 
and tx,t2,. • • ,tfc are destinations. H is the forward channel matrix from 
sources to destinations, U is the channel matrix among sources, and W is 
the channel matrix among destinations. 



Ui[n] and fi[n] are the signal received by destination i, ; and 
source Si at time n, Zi[n] and Si[n] are the channel noise, 
Hij is channel coefficient from Sj to ti and also the channel 
coefficient from ti to Sj due to channel reciprocity, Wij is 
the channel coefficient from rj to Ti, and Uij is the channel 
coefficient from Sj to s,. 

III. Motivating Examples 

The main theme of this paper is to show how interaction 
among destinations and sources can help do interference 
alignment/neutralization and thus make concurrent wireless 
communication feasible. 

In this section, we start with a simple example of 3- 
user half duplex interference channel where communication is 
impossible without interaction from destinations, and describe 
a simple three-phase transmission scheme which makes con- 
current communication feasible by exploiting the interaction 
from destinations to sources. Then, we generalize this three- 
phase transmission scheme to a class of If -user interference 
channels where the forward channel matrix H can be written 
as the sum of a diagonal matrix and a rank 1 matrix. 

A. A symmetric 3-user interference channel 

Consider a 3-user interference channel with forward channel 
matrix 




(6) 



where all interference link gains are 1, and all direct link gains 
are 0. 

Since the direct links have gains 0, it is trivial to see that 
without interaction from destinations, no reliable communi- 
cation is possible. However, if we allow destinations to talk 
back to sources using the reciprocal channel, we can achieve 
| degrees of freedom as seen below. 

Consider the following three-phase interactive transmission 
scheme: 



yi = x 2 + x 3 + noise, 
2/2 = %i + 2:3 + noise, 
y 3 = xx + x 2 + noise. 



(7) 
(8) 
(9) 



• Phase 2: After receiving signal 2/1,2/2, 2/3 from sources in 
Phase 1, all destinations simultaneously send out — to 
sources using the reciprocal channel, and sources get 

/1 = -J/2 - 2/3 + noise 

= — 2xi — x 2 — X3 + noise, (10) 
f 2 = -yi - y 3 + noise 

= — 2x 2 — xi — 23 + noise, (11) 
h = -yi - 2/2 + noise 

= —2x3 — xi — x 2 + noise. (12) 

• Phase 3: Now each source has two sets of signals x.^ 
and fi, and each Si can simultaneously send out (— 3x^ — 
fi) to destinations using the forward channel, destinations 
will get 

y[ = (-3x 2 - f 2 ) + (-3x 3 - / 3 ) + noise 

= 2xi + noise, (13) 

y 2 = (-3xi - /1) + (-3x 3 - / 3 ) + noise 

= 2x2 + noise, (14) 

2/3 = (-3 x i - A) + {-3 X 2 - $2) + noise 

= 2x3 + noise. (15) 

Therefore, after two forward transmissions and one reverse 
transmission, each destination ti will get the desired symbol 
Xi from the corresponding source Sj without any interference. 
If we do not count the reverse transmission, the total degrees 
of freedom achieved are |. 

B. A symmetric K-user interference channel with special 
channel matrix 

We can extend the above result to if -user interference 
channel where the channel matrix is an all-one matrix except 
all diagonal entries being zero for any K > 3. 

Consider a if -user interference channel where the channel 
matrix H satisfies 



H zi = 0,V1 <i< K. 



(16) 
(17) 



We can show that the following three-phase interactive 
transmission scheme with two forward transmissions and one 
reverse transmission can make each destination ti get the 
desired symbol from source for all 1 < i < K . 

• Phase 1: All sources send their independent symbols 
xi, X2, . . . , xk simultaneously. Then each destination ti 
get 



(18) 



Phase 2: After receiving signal yi , y 2 , ■ ■ ■ , Uk from 
sources in Phase 1, all destinations simultaneously send 
out — yi to sources using the reciprocal channel. Then 
each source si gets 
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fi = 



yj + noise 

l<j<K,j& 



-(K-l) Xi - YI (K-2) Xj 



(19) 



Phase 3: Now each source Si has two sets of signals Xi 
and fi, and each Si can simultaneously send out 



(-(K 2 - 3K + 3) Xi - f t ) 



(20) 



to destinations using the forward channel. Then each 
destination ti will get 



noise 



y'i = E (-(^ 2 - 3K + 3 )*< - /*) 

l<j<K,jjti 

= (K - 2)(K - l)xi+ noise. (21) 

For K > 3, (K-2)(K-1) 7^ 0. Therefore, each destination 
^ will get the desired symbol from the source Si without any 
interference for all 1 < i < K. 

In fact, we can show that if H G (£ KxK can be written as 
the sum of diagonal matrix and a rank 1 matrix, then in general 
two forward transmissions and one reverse transmission can 
make each destination i; get the desired symbol from source 
s % for all 1 < i < K . 

Theorem 1. Consider a special K-user interference channel 
with reciprocal feedback channel (G = H T ), where the 
forward channel matrix H can be written as the sum of a 
diagonal matrix and a rank 1 matrix, i.e., 



H = D + uv , 



(22) 



where D G (£ KxK is a diagonal matrix and u, v G C Kxl are 
column vectors. 

If D is invertible, each component of u and v is nonzero, 
and 



a = 3 + 3v T D- 1 u + {v r D^ L uf ^ 1, 



Tn-l. .\2 



(23) 



then there exists a three-phase interactive transmission scheme 
such that two forward transmissions and one reverse transmis- 
sion can make each destination ti get the desired symbol from 
source Si without any interference for all 1 < i < K. 

Proof: See Appendix |A] ■ 
IV. Communication Scheme and Interference 

ALIGNMENT CONDITIONS 

The three-phase interactive transmission scheme described 
in Section [ill] is a natural scheme to exploit the interaction 
from destinations to sources to make concurrent wireless 
communication feasible. We can generalize this transmission 
scheme in the following: 





Fig. 3. Phase 1: All sources send their independent symbols simultaneously. 
And destinations get y = Hx + noise. 



• Phase 1 (forward transmission): All sources send their 
independent symbols simultaneously. And destinations 
get y — Hx + n, where n is the additive noise of the 
channel. 

• Phase 2 (interaction from destinations): After receiving 
signals from sources in phase 1, all destinations scale 
y and send back to sources using the reverse channel, 
sources get / = GDiy + n, where n is the additive noise 
of the channel. Since each source and each destination 
only knows what signals they sent out and received, the 
coding matrix D\ has to be diagonal. 

• Phase 3 (forward transmission): Now each source has two 
sets of signals x and /, and sources can send out a linear 
combination of x and / to destinations via the forward 
channel, destinations get 

y' = H(D 2 x + D 3 f) + h (24) 
= (HD 2 + HD 3 GDiH)x + HD 3 GD in + HD 3 fi + h. 

(25) 

Again, the matrices D 2 , D 3 are constrained to be diago- 
nal. 

In the above scheme, since each source and each destination 
only know what signals they sent out and received, the coding 
matrices Di,D 2 and D 3 have to be diagonal, i.e., each node 
can only do coding over what signals it sent out and received. 

To make destinations be able to decode desired signals, one 
possibility is that we choose diagonal matrices D\ , D 2 and D 3 
such that (HD 2 +HD 3 GDiH) is a diagonal matrix and every 
diagonal entry is nonzero, which makes each destination get 
the desired signal from y' without any interference. In fact we 
have used coding matrices with this property for the special 
class of interference channels presented in Section [Ell] So a 
natural question to ask is whether there always exists such 
D U D 2 and D 3 . 

Requiring that all off-diagonal entries be zero gives rise 
to K(K — 1) polynomial equations, and each diagonal entry 
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Fig. 4. Phase 2: After receiving signals from sources in phasel, all 
destinations scale y and send back to sources using the feedback channel. 
Sources get / = GD±y + noise. 
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y' = H(D 2 + D 3 f) 



Fig. 5. Phase 3: Now each source has two sets of signals x and /, and 
sources can send out a linear combination of x and / to destinations via the 
forward channel. Destinations get y 1 = H(D2X + Dgf) + noise. 



being nonzero leads to K inequalities, whereas in total there 
are only 3K variables, corresponding to the K variables in 
each of the diagonal coding matrices £>i,£>2 and D3. One 
may therefore be led to conjecture that if the number of 
variables is more than the number of equations, a solution 
should exist that simultaneously satisfies all equations and 
inqualities. Unfortunately, this conjecture is not true: there are 
cases where even though the number of variables is more than 
equations, due to the coupled constraints on the equations and 
inequalities, all off-diagonal entries of (HD2 + HD3GD1H) 
being zero implies that at least one diagonal entry is zero. 
A concrete example is the 3-user interference channel with 
reciprocal reverse channel, which will be discussed in Section 



One way to overcome this problem is to let the receivers 
decode the signal based on both y (which the destinations 
received during the first phase) and y' instead of relying 
purely on y' alone (as done in the above strategy). Now, the 
destinations can cancel out all interference in y' by taking a 
linear combination of y and y' if the interference terms in y' 
are aligned with the interference terms in y. 

More precisely, let B = HD 2 + HD 3 GDiH. Then for 
destination L, in Phase 1 it receives 



yi = HuXi + HijXj + noise, 
and in Phase 3 it receives 



y, 



If for all 1 < i < K, we have 



Hi 



or equivalently 



Bi. 



'!; r ■ vj.jv*, 



= h, Vj ^ i, 



(26) 



(27) 



(28) 



(29) 



for some constant Aj, then interference terms of y and y' for 
each destination are aligned with each other. So destination ti 
can compute y t — Ajj4 and get 



Therefore, if (Hi, 



(Ha - \BiijXi + noise. 
: — XiBu) ^ 0, or equivalently, 

Ba , Bi 



Hi 



Hi 



Vj + i, 



(30) 



(31) 



then destination ti can get the desired signal sent by source 
Si without any interference. 

Without introducing the auxiliary variables Ai, aligning 
(K — 1) interference terms for each destination corresponds 
to (K — 2) equations, and thus in total we have K(K — 2) in- 
terference alignment equations. Preserving the desired signals 
after canceling all interference terms leads to K inequalities. 
Therefore, to make the three-phase communication scheme 
work and thus achieve the optimal ^ degrees of freedom, 
we need to solve K(K — 2) equations with 3K variables and 
check that whether the solution satisfies the K inequalities. 

In some cases, we can reduce the polynomial equations 
corresponding to interference alignment equality constraints 
to linear equations, the solution of which has a closed-form 
expression and thus it can be easily verified whether the 
inequality constraints are satisfied. However, in general the 
system of polynomial equations is nonlinear, and it may not 
have closed-form solution, making it hard to check whether the 
inequality constraint can be satisfied. Our main approach is to 
convert the system of polynomial equations and polynomial 
inequalities to a system of polynomial equations, and then 
use tools from algebraic geometry to check the existence of 
solution for the system of polynomial equations. 



First, we show how to reduce the problem to checking the 
existence of solutions to a system of polynomial equations. 
Then we will present our main technique for solving the 
polynomial system using tools from algebraic geometry in 
Section [Vj 

Suppose there are N polynomial equation constraints fi = 
0, 1 < i < N, and M polynomial inequality constraints 
Qj ^ 0, 1 < j < M, over S variables {d%, di, ■ ■ ■ , ds}- By 
introducing an auxiliary variable t, we define a polynomial 
function g as 



M 



(32) 



Lemma 1. There exists (d\, di, ■ ■ ■ , ds) satisfying equality 
and inequality constraints 

fi = 0, VI < i < N (33) 

9j¥=0, Vl<i<M, (34) 

if and only if there exists (dx, di, ■ ■ ■ , ds, t) satisfying 

fi = 0, VI < i < N (35) 

5 = 0. (36) 



Proof: If (dx,di, . . . , ds) is a solution to < |33] > and ( |34| l, 
then (dx, di, . . . ,ds,t) is a solution to ( [35] l and ((36), where 

*= ( 37 > 

On the other hand, if (dx, di, ■ ■ ■ ,ds,t) is a solution to i 
and ( [361 ), men (<^i? di,.. . , ^s) satisfies both ( |33| l and ( |3"4"j ). 

Therefore, due to Lemma [TJ the problem of checking the 
existence of solutions to a system of polynomial equations 
and inequalities can be reduced to the problem of checking 
existence of solutions of a system of polynomial equations, 
which is well studied for algebraically closed field in algebraic 
geometry ifTTI . 

V. General Solution Methodology 

In this section we present our main technical results on 
checking the existence of solutions to the interference align- 
ment equations for generic channel parameters H and G using 
tools in algebraic geometry. For more details on algebraic 
geometry, we refer the reader to Appendix [B] and the excellent 
textbook IfTTI . 

We follow the standard notation in IfTTI . Let k denote a field, 
and let fc[£i, • • • i £n] denote the set of all polynomials in 
the variables £1, £2, ■ ■ • , £n with coefficients in k. 

Definition 1 (Definition 1 on page 5 of IfTTI ). Let k be afield, 
and let fx, . . . , f s be polynomials in fc[£i, £2, ■ • • , Then we 
call 

V(fi, ...,f s ) = {(ai, ...,a n )ek n I /;(ai, . . .,a n ) = 0, 

VI < i < s} (38) 

the affine variety defined by fx, . . . , f s . 



As discussed in Section IV our problem can be reduced to 
checking the existence of solutions to a system of polynomial 
equations. In the language of algebraic geometry, the prob- 
lem is to check whether the affine variety defined by some 
polynomials is an empty set or not, which is well studied for 
algebraically closed field in algebraic geometry. In wireless 
communication, the channel coefficients are represented as 
complex numbers C, which is an algebraically closed field. 

The standard approach to checking whether an affine variety 
is an empty set is to use Buchberger's algorithm to compute the 
Grobner basis of the given polynomials, and from the Grobner 
basis we can easily conclude whether the corresponding affine 
variety is empty or not tT2l . 

One important implication of these results in algebraic 
geometry is that if the coefficients in the polynomial equations 
are rational functions of variables {hi}, then except a set 
of {hi} which satisfies a nontrivial polynomial equation on 
{hi} and thus has a measure 0, either for all numeric values 
of {hi} the polynomial equations have a solution, or for all 
numeric values of {hi} the polynomial equations do not have 
a solution. 

More precisely, 

Theorem 2. Let fx,...,f s be polynomials in fc[£i, £2, ■ • ■ , £n], 
where k = C and all coefficients of the polynomials are 
rational functions of variables hx, h,2, • ■ • , h m . Then there 
exists a nontrivial polynomial equation on hx, h%, . . . , h m , 
denoted by e(hx, hi, ■ ■ • , h m ), such that except the set of 
{(hx, /i2> • • • > h m ) j e(hx, hi, . . . , h m ) = 0}, either for all 



(hi, h^, 
(hi, h 2 , 



■ ■ , h m ) € 



e C lxm , V(fi, 
lxm , V(fi,.. 



■■Js) * 

,/ s ) = 0. 



or for all 



Proof: See Appendix [C] ■ 
For the polynomial equations describing the interference 
alignment problem, the coefficients of the polynomials are 
rational functions of channel parameters H and G in symbolic 
form. Therefore, in the context of the interference alignment 
feasibility problem, this main result can be restated as follows. 
Either one of the following two statements hold: 

• For almost al[] channel realizations of H and G, there 
exists solution to the system of interference alignment 
equations. 

• Or, for almost all channel realizations of H and G, 
there does not exist solution to the system of interference 
alignment equations. 

Although in theory we can compute in finite number of 
steps the symbolic Grobner basis for these polynomials with 
symbolic coefficients to check whether for almost all H and 
G there exists solutions, it turns out to be computationally 
infeasible to run the Buchberger's algorithm for the symbolic 
polynomial equations for most of our interference alignment 
problems, due to the fact that the orders of intermediate 
symbolic coefficients can increase exponentially. However, it 

'We emphasize that in our results "almost all" means for all numeric values 
except a set of parameters which satisfy a nontrivial polynomial equation. 
Therefore, in contrast to results in |5 ] and [4], our results are not sensitive to 
whether the channel parameters are rational or irrational. 



is much easier to compute a Grobner basis numerically. Due 
to Theorem [2] 

Corollary 1. If we draw the channel parameters according to 
a continuous probability distribution, then with probability one 
the numeric polynomial equations have a solution if and only 
if for almost all channel realizations the polynomials equations 
have solution. 

Hence, while we may not be able to prove that certain 
polynomial equations have solution for almost all channel 
parameters because of computational difficulty, numeric sim- 
ulations can let us make claims with high credibility. 

VI. Interactive Interference Alignment for 

SMALL K WITH HALF DUPLEX RADIOS 

In this section, we study interactive interference alignment 
for half-duplex if -user interference channel with small K. 

A. K = 3 



As discussed in Section IV to make the three-phase in- 
teractive communication scheme work, ideally we would like 
to find diagonal coding matrices D\, D 2 and D 3 to make 
(HD2+HD3GD1H) be a diagonal matrix with each diagonal 
entry being nonzero, so that interference neutralization can be 
achieved in the third phase. For 3-user interference channel, 
the number of interference neutralization equations is 6, and 
the number of variables is 3K = 9. While it seems we 
can solve the interference neutralization equations with more 
variables, due to the coupled constraints on the equations and 
inequalities, all off-diagonal entries of (HD2 + HD3GD1H) 
being zero implies that at least one diagonal entry is zero. A 
simple counter example is the case when the reverse channel 
is reciprocal to the forward channel, i.e., G = H T . We did 
extensive numeric simulations, and in each numeric instance 
we find that there is no solution to solve both equations and 
inequalities by using the methodology (computing Grobner 
basis) introduced in section [V] Therefore, due to Corollary [Tj 
we believe the following claim holds: 

Claim 1. For 3-user interference channel with reciprocal 
feedback channel (G = H T ), for generic matri^ H, there 
does not exist Dx, D2 and D3 such that (HD2+HD3GD1H) 
is a diagonal matrix and every diagonal entry is nonzero. 

As discussed in Section |IV| it is not necessary to do exact 
interference neutralization. Instead, interference alignment is 
also sufficient. When K = 3, the number of variables is 
3K = 9, which is much more than the number of interference 
alignment equations K(K — 2) = 3, so it is not surprising that 
we can find D\ , D2 and D3 to solve the interference alignment 
equations and inequalities. [ 1 1 gives an analytical solution to 
do interference alignment for 3-user interference channel with 
feedback. Here, we present a different proof technique, which 
also will be applied in interference alignment with full-duplex 
antennas in Section fVTll 



First we state a lemma which will be used in the proof of 
Theorem [3] 

Lemma 2. Let So C C™ be a linear subspace with dimension 
c?q. For any m linear subspaces Si, S 2 , . . . , S m C So with 
dimensions d\ , . . . , d m respectively, if 



di < do, VI < i < m, 
then So\{U^L 1 Si) is a nonempty set. 



(39) 



Proof: Let ei, . . . ,ed G C™ be a set of bases of Sq. 
Let A e C d ° xd ° m such that every d by d submatrix of A 
has full rank. It is easy to show the existence of A with such 
property. Indeed, if every entry of A is independently gener- 
ated according to Gaussian distribution (or other continuous 
probability distributions), then with probability 1 every do by 
do submatrix of A has full rank. 

Define d Q m vectors bi,b 2 , ■ ■ ■ , bd om by setting 



3=1 



(40) 



for all 1 < i < dom. Since every do by do square submatrix 
of A has full rank, every do vectors of b\, 62, • • ■ , bd a m are 
linearly independent. 

Suppose S'o\(U"l is a nonempty set, i.e., Sq — U'^Si. 
Then b\, b 2 , 6d m € U™^. By pigeonhole principle, 
there exists i* < dom such that Si* contains at least do vectors 
of bi,b 2 , b dom . Since every d vectors of 61, b 2 , ■ ■ ■ , b dom 
are linearly independent, the dimension of Si* is no less than 
do, contradicting with the fact that di» < do- 

Therefore, 5 \(U™ 1 S i ) is a nonempty set. 

■ 

Theorem 3. For 3-user interference channel with reciprocal 
feedback channel (G = H T ), and for generic channel matrix 
H, there exists solutions to the interference alignment equa- 
tions and inequalities. 

Proof: The K(K — 2) = 3 equations are 



(41) 
(42) 
(43) 

(44) 
(45) 
(46) 



B12 


_ -B13 


H\2 


H13 


B21 


B23 


H21 


H23 


B31 


B32 


H31 


H32 



and the K = 3 inequalities are 







B\2 


H u 


H12 


B22 




B23 


H22 


H23 


B33 




B32 


H33 


H 32 



2 Generic matrix means all numeric matrices except a set of matrices the 
entries of which satisfy a nontrivial polynomial equation. 



Since the number of variables is much more than the number 
of equations, we do not need to use all the variables. In fact, 



we can set D x to be a simple diagonal matrix 




Write Z?2 and D% as 



D* 



D21 











D22 











D23 


D 31 











D 32 











£>33 



(47) 



(48) 



(49) 



Then equations (41) , |42} and ( |4"3") are homogeneous linear 
equations in terms of variables D 2 i, D 2 2, D23, -D31, -D32 and 
-D33. Let 5 denote the set of solutions to linear equations 
( pT) , ( pE) and ( |43) . It is easy to check that the corresponding 
coefficient matrix has a rank of 3. Therefore, S is a linear 
subspace of dimension 6 — 3 = 3. 

Consider the following equations corresponding to the in- 
equality constraints 



(50) 





_ B12 


H n 


H\2 


B22 


B23 


H22 


H23 


B33 


B32 


H33 


a 32 



(51) 
(52) 



Let S-y denote the set of solutions to (j4TJ, ((42}, (|43} and 
( |5()) , let S 2 denote the set of solutions to (|4T}, |42}, (43} and 
and let S3 denote the set of solutions to (|4T}, |42}, ( |43) 
We can check that all the coefficient matrices of the 



(51 



and 

linear equations have a rank of 4, so all of S\ , S2 and S3 are 
linear subspaces of dimension 6 — 4 = 2. 

Now consider the set S* = S\(Si US 2 US 3 ). Any element 
in S* satisfies both the equations ( |4T) , |42}, |43}, and the 
inequality constraints ( |44) , ( |4"5) , ( ftp") . 

Since 5 is linear subspace of dimension 3, which is larger 
than the dimensions of Si, S2 and S3, S* is a nonempty set 
by Lemma [2] This completes the proof. 

■ 

The same result continues to hold when the feedback 
channel matrix is also a generic matrix independent of if. 

It is worth pointing out that the 3 -user interference channel 
under an alternate local feedback model, where each receiver 
sends feedback to its corresponding transmitter, was studied 
in lfT~3l recently. In contrast to our model, where obtaining 
feedback using natural reciprocal channel takes up a single 
time slot, local feedback can take 3 time slots to implement. 

B. K = A 

For 4-user interference channel, we have K(K — 2) = 8 
polynomial equation constraints and K = 4 inequality con- 
straints with 3K = 12 variables. We cannot use the same 
technique as in the case of three-user interference channel by 
setting D\ or D 3 to be a deterministic constant matrix and 



thus reducing the polynomial equations to linear equations. 
The reason is that if we do not use Di or D 3 as variables, 
the number of equations will be the same as the number of 
variables, and this will prevent us from finding a solution 
which also satisfies the inequality constraints. Instead, we use 
the tools introduced in Section [V] to study the feasibility of 
interactive interference alignment. More precisely, we first 
convert the system of polynomial equations and polynomial 
inequalities to a system of polynomial equations, and then 
we check whether it has solutions or note by computing the 
Grobner basis of corresponding polynomials. 

Due to computational difficulty, it is hard to compute the 
Grobner basis for the symbolic polynomials. We do extensive 
numeric verifications: each time we generate random numeric 
channel matrices, and then apply Buchberger's algorithm to 
compute a Grobner basis. In all numeric instances, the Grobner 
basis does not contain a nonzero scalar, which implies that 
solution to ((28} and ( (3T) exists. Therefore, due to Corollary 
[T] although we cannot compute a symbolic Grobner basis 
explicitly, from numeric verifications we believe the following 
claims hold. 

Claim 2. For 4-user interference channel with reciprocal feed- 
back channel (G = H T ), for generic matrix if, there exists 
Di, D 2 and D 3 such that each row of {HD 2 + HD3GD1H) 
is proportional to the corresponding row of if except the 
diagonal entries. 

Claim 3. For 4-user interference channel with out-of-band 
feedback channel, for generic matrices if and G, there exists 
D\, £>2 ond D3 such that each row of (HD2 + HD 3 GD\H) 
is proportional to the corresponding row of H except the 
diagonal entries. 

[ 10 1 formulates the interference alignment problem as a rank 
constrained nonconvex optimization problem, and reports that 
all numeric simulations confirm the existence of solutions to 
the interference alignment equations and inequalities. 

In some cases when H has special structures, we are able 
to compute the symbolic Grobner basis of the interference 
alignment polynomials, and thus we can conclude whether 
for almost all H with such structure there exists interference 
alignment solutions. One positive example is that if if is a 
symmetric matrix and all diagonal entries are zero, then for 
almost all such if, the interference alignment equations are 
solvable. 

Theorem 4. For 4-user interference channel with reciprocal 
feedback channel (G = H T ), if if is symmetric channel and 
all diagonal entries are zero, i.e., if can be written as 



( 



if = 



V 




hi 

h 3 



hi 


h 5 



h 2 
hi 


/'6 



h 3 
hn 
h e 




/ 



(53) 



then for generic {hi,h2,h3,h4,h^,hQ}, i.e., for all 
{hi, /12, /13, hi, /15, he} except a set of {hi, hi,, h%, hi, /15, he} 
which satisfies a nontrivial polynomial equation, there 



exists Di,D 2 an d D 3 such that (HD2 + HD 3 GDiH) 
is proportional to the corresponding row of H except the 
diagonal entries. 

Proof: We derive the interference alignment equa- 
tions and then use commercial computer algebra system 
Maple to compute the corresponding symbolic Grobner basis, 
which does not contain a nonzero scalar. Therefore, due to 
Corollary [4] in Appendix [B] we conclude that for generic 
{hi, fi2, h 3 , li4, h$, he}, there exists D\, D2 and D 3 such that 
(HD2 + HD3GD1H) is proportional to the corresponding 
row of H except the diagonal entries. ■ 

C. K — 5 and K = 6 

Recall that the number of variables N v is 3K, the number of 
equations N e is K(K — 2) and the number of inequalities iVj e 
is K. N e grows faster than N v as K increases, and when K > 
5, N v < N e . We expect that when the number of variables is 
not more than the number of inequations, a valid solution may 
not exist. Indeed, we confirm this conjecture for K — 5 and 
6 via computing Grobner basis for random numeric H and G 
and finding that 1 is in the Grobner basis in all instances. 

Claim 4. For K-user interference channel with reciprocal 
feedback channel (G — H T ) where K = 5 or 6, for generic 
matrix H, there does not exist D± , D2 and D 3 such that 
each row of (HD 2 + H D 3 GDiH) is proportional to the 
corresponding row of H except the diagonal entries. 

Claim 5. For K-user interference channel with out-of-band 
feedback feedback channel where K = 5 or 6, for generic 
matrices H and G, there does not exist D\,D2 and D 3 such 
that each row of (HD2 + HD 3 GD\H) is proportional to the 
corresponding row of H except the diagonal entries. 

Since one reverse transmission from destinations to sources 
is not enough, in Phase 2 (interaction phase) we can let 
destinations make more than one reverse transmissions to 
sources. In one extreme case, if each destination sends what 
it receives in Phase 1 to all sources sequentially, where in 
total we use K reverse transmissions^) then each source can 
compute the symbol signals sent out by all sources in the first 
phase, and then each source can do precoding over all symbols 
in the second phase, which reduces to a MIMO broadcast 
channel. 

It turns out that when K — 5 or K = 6, two reverse trans- 
missions in Phase 2 is sufficient to do interference alignment 
while still preserving the desired signals. A natural three-phase 
scheme is as follows: 

• Phase 1: All sources send their independent symbols 
simultaneously. And destinations get y = Hx + n, where 
n is the additive noise of the channel. 

• Phase 2-1: After receiving signals from sources in phase 1, 
all destinations scale y and send back to sources using 
the feedback channel, sources get f 1 = GD\y + noise. 

3 In fact, it is easy to see that (K — 1) reverse transmissions are also 
sufficient. 



• Phase 2-2: All destinations scale y using a different 
coding matrix and send back to sources, sources get 
f 2 = GD 2 y + noise. 

• Phase 3: Now each source has three sets of signals x, f 1 
and f 2 , and sources can send out a linear combination 
of x, f 1 and f 2 to destinations via the forward channel, 
destinations get 

y' = H(D 3 x + DJ 1 + D 5 f 2 ) + noise (54) 

= (HD 3 + HD4GD1H + HD 5 GD 2 H)x + noise. 

(55) 

With two reverse transmissions, the number of variables is 
increased to hK from 3K, and is more than N e = K(K — 2) 
for K = 5, 6. Extensive numeric simulations show that for 
generic matrix H and G, there exists diagonal coding matrices 
{-Dj}i<i<5 to align interferences while still preserving desired 
signals. 

Claim 6. For K-user interference channel with reciprocal 
feedback channel (G — H T ) where K = 5 or 6, for 
generic matrices H, there exists Di,D2,D 3 ,D^ and -D5 
such that each row of (HD 3 + HD 4 GD 1 H + HD 5 GD 2 H) 
is proportional to the corresponding row of H except the 
diagonal entries. 

Claim 7. For K-user interference channel with out-of-band 
feedback channel where K = 5 or 6, for generic matrices H 
and G, there exists D\, D2, D 3 , D 4 and D§ such that each 
row of (HD 3 + HDiGDxH + HD 5 GD 2 H) is proportional 
to the corresponding row of H except the diagonal entries. 

VII. Interactive Interference Alignment for 

SMALL K WITH FULL DUPLEX RADIOS 

In this section, we study how to exploit interaction in 
interference channel where all sources and destinations have 
full-duplex antennas so that all nodes can transmit and receive 
in the same band simultaneously. We will show how full- 
duplex operations can help align interferences and achieve the 
optimal degrees of freedom for some classes of interference 
channels. 

A. SISO IFC with small K 

Recall in the system model of full-duplex IFC, H is the 
channel matrix from sources to destination, U is the channel 
matrix among sources, and W is the channel matrix among 
destinations. Consider the following simple two-phase trans- 
mission scheme. 

• Phase 1: All sources send out signals x simultaneously. 
After the transmission, sources get / = Ux + noise, and 
destinations get y = Hx + noise. 

• Phase 2: sources send out a linear combination of x and /, 
and destinations send out scaled version of y. Therefore, 
destinations get H(D\x + D 2 f) + WD 3 y + noise = 
(HD 1 +HD 2 U + WD 3 H)x + noise, where D U D 2 and 
D 3 are diagonal coding matrices. 

Note that in this transmission scheme no feedback channel 
from destinations to sources is required. Similarly, if each 



row of (HD X + HD 2 U + WD 3 H) is proportional to the 
corresponding row of H except all diagonal entries, then 
interferences at all destinations are aligned and all destinations 
can retrieve their designed signals from sources without any 
interference, and thus achieve the optimal y degrees of 
freedom. 

Our first result is that the above two-phase transmission 
scheme works for K — 3 and K = 4. 

Theorem 5. For K — 3 and K = 4, for generic channel 
matrices U, H and W, there exists diagonal matrices D\ , D 2 
and D 3 such that each row of (HDi + HD 2 U + WD 3 H) 
is proportional to the corresponding row of H except all 
diagonal entries. 

Proof: The interference alignment equations are linear in 
Di, D 2 and D 3 , so as in the proof of Theorem [3] we can use 
dimension argument to show that there always exists a solution 
to the system of linear interference alignment equations and 
it also satisfies the inequality constraints. ■ 
Using the same dimension argument, we can prove that the 
above scheme does not work for K = 5 or bigger. 

B. MIMO IFC with K = 4 

Next we show that for four-user MIMO interference chan- 
nel, full-duplex antennas can help design a simple transmission 
scheme to achieve the optimal degrees of freedom. 

For if -user MIMO interference channel where all sources 
and destinations are equipped with M antennas, [1] shows 
that for K = 3, M > 1, vector space interference alignment 
can achieve the optimal = degrees of freedom and 
channel diversity is not required. [3] and [14] prove that in 
general vector space MIMO interference alignment can at 
most get j^pf degrees of freedom, which is strictly less 
than the optimal when K > 4. We show that for four- 
user MIMO interference channel with full-duplex antennas, a 
simple two-phase transmission scheme can achieve the optimal 
2M degrees of freedom at least for all M < 15. 

Suppose source Sj wants to send M symbols 
B(i-i)M+uX(i-i)M+2>--->ZiM t0 destination t. t , for all 
1 < i < 4. Let x = (xi,x 2 , . . . ,X4m) T - Consider a natural 
two-phase transmission scheme for 4-user MIMO IFC with 
full-duplex antennas: 

• Phase 1: All sources send out signals x simultaneously, 
i.e., the ith antenna sends xi for all 1 < % < 4M. 
After the transmission, sources get / = Ux + noise, and 
destinations get y 1 = Hx + noise. 

• Phase 2: sources do linear coding over x and /, and send 
D\x + D 2 f, and destinations do linear coding over y 1 
and send D 3 y, where D%, D 2 and D 3 are block diagonal 
matrices with diagonal block size of M x M . Therefore, 
destinations get 

y 2 = H(D lX + D 2 f) + WD 3 y x + noise 

= (HD 1 +HD 2 U + WD 3 H)x + noise. (56) 

Note that since each node can only do coding over what 
signals it has or received from other nodes, the coding matrices 



Di , D 2 and D 3 have to be block diagonal matrices with 
diagonal block size of M x M. 

A sufficient condition for all destinations to retrieve the 
desired symbols from the corresponding sources without in- 
terference is that at each antenna of destination U, all in- 
terference symbols are aligned except the desired symbols 
X{i-i)M+i,%{i-i)M+2i ■ ■ -,x iM , for all 1 < i < 4. 

More precisely, let B = HD 1 + HD 2 U + WD 3 H and 
denote the index set {(i - l)M + 1, (i - l)M + 2,..., iM} 
by Ii. Then for the destination ti, in Phase 1 it receives 

vln = ^2 H mm /x' m + ^2 HmjXj + noise term, (57) 



and in Phase 2 it receives 



B 



mm' ^rn 



noise term. 



(58) 



The first term on the RHS of ( |57| ) and ( |57| i correspond to 
the desired signals, and the second term are interferences. If 
at each antenna of ti, all interferences are aligned, i.e., for all 
m € Ii, there exists A m G C, such that 



H„ 



rnj 



then by computing y n 
interferences 



Y? i It, (59) 
- Xmy^, U can cancel all 



Vn 



^rnVm 



noise term. 



(60) 
(61) 



Now ti gets k linear equations in k variables {x m } me j i . If 
the coefficient matrix of the k linear equations has full rank, 
then ti can retrieve all desired signals {x m } m ^i i corrupted by 
noise only. 

The condition that the coefficient matrix has full rank is 
equivalent to the condition that the determinant of the coeffi- 
cient matrix is nonzero. Therefore, to make the transmission 
scheme work, we need to find block diagonal coding matrices 
Di, D 2 , D 3 to satisfy the interference alignment linear equa- 
tions ( |59"l l and make the determinants of coefficient matrices be 
nonzero. Without introducing the auxiliary variables A m , we 
have KM((K — 1)M — 1) interference alignment equations, 
which are linear in terms of Di,D 2 and D 3 . Since the 
equations are linear, there exists analytical solution to D\ , D 2 
and D 3 which are rational functions of U, H and W. Then, 
we can plug the analytical solution into the expression of the 
determinant of the coefficient matrix and check whether it is 
zero. We have checked that for K = 4, the determinants of the 
coefficient matrices are nonzero expressions for all M < 15. 

Theorem 6. For 4-user MIMO interference channel where 
each node have M full-duplex antennas, then for generic 
matrices H, U and W, the two-phase scheme can achieve the 
optimal y degrees of freedom, at least for M < 15. 



3-User Interference Channel 




Fig. 6. Numeric Finite SNR Performance Simulations for 3-User Interference 
Channel with Out-of-Band Reverse (Feedback) Channel. 



4-User Interference Channel 




Fig. 7. Numeric Finite SNR Performance Simulations for 4-User Interference 
Channel with Out-of-Band Reverse (Feedback) Channel. 



Note that for two-phase scheme, in total we have N v — 
2>KM 2 variables and N e = KM((K - l)M - 1) linear equa- 
tions. When K = 4, N v = 12M 2 > N e = 12M 2 - 4M for 
all M, which means we have more variables than equations, 
and thus we believe Theorem [6] holds for any M. 

VIII. Numeric Finite SNR Performance Evaluation 

In this section, we present simulation results on the finite 
SNR performance of our three-phase interactive interference 
alignment scheme with out-of-band reverse (feedback) trans- 
missions for 3-user and 4-user interference channels, where 
the channel parameters are drawn from i.i.d. complex Gaussian 
distribution. We assume both sources and destinations have the 
same power constraint. We formulate the problem as a non- 
convex optimization problem, and use generic local algorithm 
to maximize the sum rates among all possible interference 
alignment solutions while satisfying power constraints. For 
each SNR, we simulate over 50 channel realizations and 
compute the average of maximum sum rates. We compare the 
maximum rates achieved via interactive interference alignment 
with the standard time sharing scheme. 

For K = 3 and K = 4, our interactive interference 
alignment scheme can achieve the optimal y degrees of 
freedom, while time sharing scheme only gets 1 degree of 
freedom. Therefore, at high SNR, we expect our scheme has 50 
percent improvement for K = 3 and 100 percent improvement 



for K — 4 over time sharing scheme, and this matches our 
simulation results in Figure [6] and Figure [7] 

IX. Discussions and Open Problems 

In this section, we give a general multi-phase interactive 
communication scheme for X-user interference channel with 
large K and discuss how 0(\/K) degrees of freedom might 
be achieved with reverse transmissions taken into account. We 
also discuss some open problems following this line of work. 

A. Multi-Phase Transmission Scheme for General K 



In Section VI we have shown that for general K we can 



achieve the optimal Ir degrees of freedom using the three- 
phase scheme with (K — 1) reverse transmissions in Phase 2. 
However, if we take reverse transmissions into account, the 
total number of transmissions is 2 + (K — 1) = K + 1, and 
thus the total degrees of freedom is j^r^ < 1. In this section, 
we give a multi-phase transmission scheme and use equation- 
variable counting argument to show that we may potentially 
achieve 0(yK) degrees of freedom even if we take reverse 
transmissions into account. 
Multi-Phase Transmission Scheme 

• Forward Phase 1: All TXs send their independent sym- 
bols simultaneously. And RXs get y 1 = Hx + noise. 

• Reverse Phase 1: All RXs scale y and send back to TXs 
using the feedback channel. TXs get f 1 = GZ?ij/+noise. 

• Forward Phase 2: All TXs send a linear combination of 
x and f 1 to RXs. And RXs get y 2 = H{D 2 x + D 3 f 1 ) + 
noise. 

• Reverse Phase 2: All RXs send a linear combination of y 1 
and y 2 to RXs. TXs get f 2 = GiD^y 1 + D 5 y 2 ) + noise. 

• ... 

• Forward Phase i: All TXs send a linear combination 
of x, f , / 2 , . . . , f 1 to RXs, using i diagonal coding 
matrices, and RXs get y % . 

• Reverse Phase i: All RXs send a linear combination of 
y x ,y 2 , . . . ,y l to RXs, using i diagonal coding matrices, 
and RXs get y l . 

• ... 

• Forward Phase N: All TXs send a linear combination 
of x, f 1 , f 2 , . . . , f 1 ^^ 1 to RXs, using N diagonal coding 
matrices, and RXs get y N . 

The number of coding matrices is 



l + 2 + 2 + 3 + 3 + -- - + (JV-l) + (JV-l)+JV 



N 2 



1, 



(62) 



so the total number of variables N v is (N 2 — 1)K. On the other 
hand, the total number of interference alignment equations N e 
is K(K-2). If we want to have N v > N e , then N > yjK - 1, 
and thus N — Cl(yK). If there exists a solution to the 
corresponding system of polynomial equations, then the sum 
degrees of freedom achievable is 



K 



< 



K 



2N ~ 2^/W^l 



= 0(VK) 



(63) 



Therefore, by equation-variable counting, we conjecture that 
the multi-phase interactive interference alignment may achieve 
0{VK) de grees of freedom for general if. 

B. Open Problems 

We list several interesting open problems following this line 
of work on interactive interference alignment. 

• Constructive Proof. While we have used the generic 
Grobner basis tool to show whether the interference 
alignment equations are solvable or not, it will be more 
interesting to derive an analytical solution, like what we 
have done in the case where H can be represented as the 
sum of a diagonal matrix and a rank 1 matrix. 

• Efficient Algorithm and Optimization. In theory we 
can use the Buchberger's algorithm to get an interference 
alignment solution, but the computational complexity is 
high. So one natural question is whether we can design 
a much more efficient algorithm to solve the interference 
alignment equations. In addition, as presented in Section 



VIII to maximize the sum rates using interactive interfer- 



ence alignment, we formulate a nonconvex optimization 
problem and use generic local algorithm to solve it. Can 
we design an efficient algorithm to solve the optimization 
problem so that it can be used in practice (in a distributed 

way)? 

• Half-duplex interactive IA for general if. We give 
an equation-variable counting argument to show why 
0{yK) degrees of freedom might be achievable using 
interactive interactive alignment with reverse transmis- 
sions into account. Can one rigorously prove or disprove 
it? Can we design a new transmissions scheme to get 
degrees of freedom linearly growing with K using the 
same idea of interactive interference alignment? 

• Full-duplex interactive IA for general if. We have 
shown for small if, e.g., K = 3 and 4, interactive IA 
with full-duplex radios can achieve the optimal degrees 
of freedom, while for general if the current two phase 
scheme does not work for if > 5. It is interesting to study 
how to extend the two phase interactive IA to general K- 
user interference channel with full-duplex antennas. 

• Interactive IA with relays. Relays naturally fit into 
our model, for both half-duplex and full-duplex radios. 
Through interaction among relays, sources and destina- 
tions, relays can help sources do interference alignment. 
For instance, for half-duplex if -user interference channel, 
with aif 2 relays, for some constant a, in the natural three 
phase scheme with one reverse transmission from the 
destinations, relays can help align all interferences at each 
destination. Can one characterize the minimum number 
of relays needed, or design a new transmission scheme 
that achieves alignment using o(if 2 ) relays? 

X. Conclusion 

We study new channel models where interaction among 
sources and destinations is enabled, e.g., both sources and 
destinations can talk to each other. The interaction can come 



in two ways: 1) for half-duplex radios, destinations can talk 
back to sources using simultaneous out-of-band transmission 
or in band half-duplex transmission; 2) for full-duplex radios, 
both sources and destinations can transmit and listen in 
the same channel simultaneously. Although [9| shows that 
for interference channel, relays, feedback, and full-duplex 
operation cannot improve the degrees of freedom beyond 
Y, we demonstrate that the interactions among sources and 
destinations enables flexibility in designing simple interference 
alignment scheme and in several cases achieves the optimal 
degrees of freedom. 

We present a three-phase interactive communication scheme 
for half-duplex if -user interference channel with small if, 
a two-phase interference communication scheme for full- 
duplex if -user interference channel with small if, and show 
that the interactive interactive alignment scheme can achieve 
the optimal degrees of freedom for several classes of IFC, 
including half-duplex 3-user IFC, full-duplex 3-user and 4- 
user IFC, and full-duplex 4-user MIMO IFC. We do extensive 
numeric simulations and show the proposed interactive inter- 
active alignment scheme also works for some other classes 
of IFC empirically. We use tools from algebraic geometry to 
show why success of numeric simulations can suggest the 
scheme should work well for almost all channel parameters 
in a rigorous way. Simulation results on the finite SNR 
performance of interactive interactive alignment for 3-user and 
4-user half-duplex interference channels are presented. Finally, 
we discuss a general multi-phase interactive communication 
scheme for if -user interference channel with large if, and list 
several open problems. 

Appendix A 
Proof of Theorem[T] 

Proof of Theorem |7J Following the notation introduced 
in Section IV equivalently, we need to prove that there exists 



diagonal coding matrices Di,D 2 and D3 such that (HD 2 + 
HD3GD1H) is a diagonal matrix and each diagonal entry is 
strictly nonzero. 

We give a constructive proof. First, we expand (HD 2 + 
HD 3 GD 1 H) and get 

HD 2 + HD 3 GD 1 H (64) 

=(D + uv T )D 2 + (D + uv T )D 3 (D + vu T )D 1 (D + uv T ) 

(65) 

=DD 2 + uv T D 2 + 

(DD 3 + uv T D 3 )(DD 1 + vu T Di)(D + uv T ) (66) 
=DD 2 + uv T D 2 + (DD 3 DD 1 + uv T D 3 DD! + 

uv T D 3 vu T D 1 + L>£> 3 uu T L>i)(T> + uv T ) (67) 
=DD 2 + uv T D 2 + DD 3 DDiD + uv T D 3 DD 1 D+ 

uv T D 3 vu T D 1 D + DD 3 vu T D 1 D+ 

DD 3 DD lU v T + uv T D 3 DD lU v T + 

uv T D 3 vu T D lU v T + DD 3 vu T D lU v T . (68) 



Define U = diag(u), i.e., U is a diagonal matrix where the 
ith diagonal entry is the ith component of u. Similarly, define 
V = diag(u). Set 

Di = D^U^V, (69) 
D 3 = D- 1 V- 1 U. (70) 

Then 

D 3 DDiD = I, (71) 
u T D\D = v T , (72) 
DD 3 v = u, (73) 

where I is the identity matrix. 
Therefore, 

HD 2 + HD 3 GDiH (74) 
=DD 2 + uv T D 2 + D + uv T + uv T D 3 vv T 

rr\ rj-f rri rj-\ 

+ uv + uv + uu D\uv + 

uv T D 3 vu T D\uv T + uu T D\uv T (75) 
=D{D 2 +I)+ uv T D 2 + (3 + 2v T D 3 v + u T D lU + 

v T D 3 vu T Diu)uv T (76) 
=D(D 2 +I)+ uv T D 2 + 

(3 + 3v T D- 1 u+(v T D- 1 u) 2 )uv T . (77) 

Define a = 3 + 3v T D~ 1 u + (w T .D -1 u) 2 , and set 

D 2 = -al. (78) 

Then 

HD 2 + HD 3 GDiH = (1 - a)D. (79) 

Since a^l, (HD 2 + HD 3 GDiH) is a diagonal matrix and 
each diagonal entry is nonzero. 

■ 

Appendix B 

Mathematical background on algebraic geometry 

In the section we present some basic terminologies and 
results in algebraic geometry. In particular, we introduce the 
Grobner basis and Buchberger's algorithm, which are standard 
tools to check whether a system of polynomial equations have 
solutions or not. For more details on algebraic geometry, we 
refer the reader to the excellent textbook ifTTI . 

Recall that we use k to denote a field, /c[£i, £2, ■ ■ • , £n] 
to denote the set of all polynomials in £1,^2, with 
coefficients in k, and use V(fi, . . . , f s ) to denote the affine 
variety defined by polynomials f t , . . . , f s G £ 2) • • • , £n]- 

Definition 2 (Definition 2 on page 30 of IfTTI ). Let fx,.. . ,f s 
be polynomials in k[^i,^ 2 , . . . ,£ n ]. Then we denote the ideal 
generated by fx, . . . , f s as 

s 

(ft, ••-,/>) - (^hift I hi,...,h s € fe[6>6> ■•• >£«]}■ 
i=l 

(80) 



It is straightforward to see that 

V(f 1 ,...,f s )=V((f 1 ,...,f s )). (81) 

Theorem 7 ( The Weak Nullstellensatz on page 170 of 
111 111 ). Let k be an algebraically closed field and let I C 
£2, • ■ ■ j £n] ^ e an ideal satisfying V(I) = 0. Then 

J = ••,£»]• 

Corollary 2. Lef fi, ■ ■ ■ , f s be polynomials in 
M£l> ^2) • • • > £«.]> where k is an algebraically closed field. 
Then V(fi, ...,f s ) = % if and only if 1 £ . . . , / s ). 

Prao/- If 1 e then (/i,...,/ s ) = 

fc[£i,...,£n], and thus F(/ x ,...,/ s ) = V((/i, . . . , /«)) = 

y(fc[£ 1 ,...,U) = 0. 

On the other hand, if V(/i, . . . , f s ) = 0, equivalently, 
V((/i, . . . , /s)) = 0, then by Theorem [7] we have 

■■■Js)= k[£i, . . .,£„], (82) 

and thus 1 e . . . ,/ s ). 

■ 

Therefore, given polynomials f u . . . , f s E k[£i, £„], 
where k is an algebraically closed field, to check whether 
V(fi, . . . , f s ) is nonempty, i.e., whether there exists solution 
to the system of polynomial equations 

/* = 0,Vl<z<s, (83) 

we only need to check whether 1 g (fi, ■ ■ ■ , fs}- This is the 
ideal membership problem, which is well studied in algebraic 
geometry. The basic tool we will use is Grobner basis |12|. 
Below we give a short introduction to Grobner basis and 
show how to use Grobner basis to solve the ideal membership 
problem. 

Definition 3 (Definition 7 on page 59 of IfTTI ). Let f = 

^2 a a a x a be a nonzero polynomial in fc[£i, . . . , £ n ] and let 
> be a monomial order. Then we define the multidegree of f 
as 

multideg(f) = max(a : a a 7^ 0), (84) 

where the maximum is taken with respect to the monomial 
order >. 

And we define the leading term of f as 

LT(f)±a muMdeg(f) r vlUde9if) - (85) 

Given an ideal /, we use LT(I) to denote the set of leading 
terms of elements of /, i.e., 

LT(I) 4 {LT{f)\f e /}. (86) 

Due to Theorem [8] every ideal in . . . , £ n ] nas a finite 
generating set. 

Theorem 8 (Hilbert Basis Theorem on page 76 of IfTTI ). Every 
ideal I C ■ • ■ 1 £n] has a finite generating set. That is, 

I = (91, ■ ■ ■ , 9t) far some gi, . . . , g t G /. 



Definition 4 (Definition 5 on page 77 of IfTTl ). Fix a monomial 
order. A finite subset G = {<?i, . . . ,g t } of an ideal I is said 
to be a Grobner basis if 

(LT( 9l ),...,LT(g t )) = (LT(I)). (87) 

Equivalently, a set {gi, . . . , gt} C I is a Grobner basis of / 
if and only if the leading term of any element of / is divisible 
by one of the LT(g t ) IfTTl . 

Proposition 1 (Proposition 1 on page 82 of IfTTl ). Let 

G = {g\, . . . , gt } be a Grobner basis for an ideal I C 
. . . , £„] anc/ fe? / £ ■ ■ • , 77ze« there is a unique 
r € fc[£i,.. . , £ n ] w/f/i the following two properties: 

(i) No term of r is divisible by any of LT(g{) 1 . . . , LT(g t ). 

( ii) There is g € I such that f = g + r. 

In particular, r is the remainder on division of f by G 
no matter how the elements of G are listed when using the 
division algorithm. 

An important property of Grobner basis is that given a 
polynomial /, we can check whether / is in the corresponding 
ideal via division of / by G. More precisely, 

Corollary 3 (Corollary 2 on page 82 of HTJ). Let G = 

{<7i, . . . , g t } be a Grobner basis for an ideal I £ • • • , £n] 
and let f £ • • ■ , £n]- Then f £ I if and only if the 

remainder on division of f by G is zero. 

Corollary 4. Let G = {gi, . . . , gt} be a Grobner basis for an 
ideal Id fc[£i, . . . , f n ]. Then 1 £ I if and only if 

9i = A, (88) 

for some i < t and nonzero A £ k. 

Proof: We first prove the easier direction. If for some 
h9i = A ^ 0, then (g t ) = k[£_ 1: . . .,£„]. So / = ...,£„] 
and thus 1 £ I. 

For the other direction, suppose that 1 £ I, but for each 
i < t, the leading term of gi has a nonzero order, i.e., gi 
does not have the form gi — A for some nonzero X £ k. 
Then after we do division of 1 by G, the remainder will be 1, 
so the remainder is nonzero, and thus 1 ^ / by Corollary [3] 
which leads to a contradiction. Therefore, for some i, gi can 
be written as gi = A for some nonzero A £ k. ■ 

Due to Corollary [2] and Corollary [4] to check whether 
there exists solution to a system of polynomial equations over 
algebraically closed field, we only need to compute a Grobner 
basis of the corresponding ideal and check whether a nonzero 
scalar is in the Grobner basis. 

The standard approach to computing a Grobner basis is 
Buchberger's algorithm (V2\ . which can be viewed as a gen- 
eralization of Gaussian elimination for linear systems. 

Theorem 9 (Buchberger's Algorithm on page 90 of ifTTIl ). 
Let I — (/i, /2, . . . , f s ) =/= {0} be a polynomial ideal. Then a 
Grobner basis for I can be constructed in a finite number of 
steps by the following algorithm: 

Buchberger's algorithm can be proven to terminate in finite 
number of steps, and one important consequence is that if the 



Algorithm 1 Buchberger's algorithm to compute Grobner 

basis 

Input: F=(f 1 J 2 ,...J s ) 

Output: a Grobner basis G — (gx,g2, ■ ■ ■ ,9t) for /, with 
F £ G 

Initialization: G £- F 
repeat 

G' <— G 

for each pair {p, q},p ^ q in G' do 

a <— multideg(p), and (3 <s— multideg(<7)) 

7 «- (71, ... , 7„), where j t = max{a % , ft) for each 

i. 

S(p, q) lt( p ) P~ LT( q ) Q ( note tnat leading terms 
will be canceled) 

q' Ql 

S <— S(p, q) , where S(p, q) the remainder 
on division of S(p, q) by G' using multivariate division 
algorithm 

if S ^ then 

G <- GUS 
end if 
end for 
until G = G' 



coefficients in the polynomial equations are rational functions 
of variables {hi}, then except a set of {hi} which satisfies 
a nontrivial polynomial equations on {hi} and thus has a 
measure 0, either for all {hi} the polynomial equations have 
a solution, or for all {hi} the polynomial equations do not 
have a solution. This result is clearly stated in Theorem [2] in 
Section fV] 

Appendix C 
Proof of Theorem|2] 

Proof of Theorem |2f Suppose all coefficients in 
fi , . . . , f s are in symbolic form, and we run Buchberger's 
algorithm to compute G = {gi, . . . , gt}, the symbolic Grobner 
basis for /a; • • • > /«)■ In eacn step of Buchberger's algo- 
rithm, we will get new derived polynomials with symbolic co- 
efficients, which are also rational functions of hi, h 2 , ■ ■ ■ , h m . 
Since Buchberger's algorithm is proven to terminate in finite 
number of steps, the set of all symbolic coefficients, which 
are all rational functions of h\,h>2, ■ ■ ■ , h m , is finite. Let e be 
the product of all numerators and denominators of the set of 
symbolic coefficients, and thus e is a nontrivial polynomial of 
hi, h,2, . . . , h m . 

For all numeric (hi, h 2 , ■ ■ ■ , h m ) £ C lxm such that 
e(hi, h 2 , ■ ■ ■ , h m ) 7^ 0, the structure of numeric Grobner basis 
is the same as the structure of symbolic Grobner basis, and 
thus either V(fi, ...,/*)= or V(fi, ...,/,) ^ 0, which 
depends on whether the symbolic Grobner basis contains a 
nonzero constant. 

Therefore, except the set of 

{(h!,h 2 , . . . ,h m ) I e(hi,h 2 ,...,h m ) = 0}, which has 
a measure 0, either for all (hi, h%, . . . , h m ) £ C lx ™\ 



V(fi,...,f s ) ? 0, or for all {h u h 2 , . . . ,h m ) e C lxm , 
V{h,...,fs)=t ■ 
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